Black Hole Instability Induced by a Magnetic Field 



Martin Ammon a , Johanna Erdmenger b , Patrick Kerner b , Migael Strydom b 

"Department of Physics and Astronomy 
University of California, Los Angeles, CA 90095, USA 
b Max- Planck- Institut fur Physik (Werner-Heisenberg-Institut) 
Fdhringer Ring 6, 80805 Miinchen, Germany 



Abstract 



In the context of gauge/gravity duality, we find a new black hole instability in 
asymptotically AdS spaces. On the field theory side, this instability is induced 
by a magnetic field in the vacuum, in contrast to previous instabilities which 
occur at finite density. On the gravity side, this corresponds to a spatial compo- 
nent of the gauge field in SU(2) Einstein- Yang-Mills theory, which provides the 
crucial non-Abelian structure. Our analysis may provide supporting evidence 
for recent QCD studies of p meson condensation induced by a magnetic field. 
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1. Introduction 

Within gauge/gravity duality, asymptotically AdS black hole solutions and 
their stability properties play a crucial role since black holes are dual to thermal 
states on the field theory side. In particular, the stability of black holes is 
connected to the stability of the thermal state. In recent years, it was found 
that many asymptotically AdS black holes are unstable, forming hair [1, 2]. 
The most prominent examples are charged AdS black holes which are unstable 
against condensation of scalar [3, 4] or gauge fields [5]. The condensate arising 
from a scalar or gauge field is a holographic description of a superfluid with 
scalar or vector condensates, respectively. 

In this letter we report on a new type of instability of Anti-de Sitter black 
holes in SU(2) Einstein- Yang-Mills theory. Our analysis is motivated by recent 
field theory studies within QCD [6, 7, 8, 9], where it has been proposed that a 
strong magnetic field may generate p meson condensation and superconductivity 
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in the QCD vacuum 1 . Essentially, in these papers it has been argued that 
in a strong magnetic field, a gluon-mediated attraction between quarks and 
antiquarks of different flavour leads to a colourless spin-triplet bound state with 
quantum numbers of an electrically charged p^ meson. These quark- antiquark 
pairs condense to form a new ground state. Within gauge/gravity duality, a 
similar result has been obtained in the Sakai-Sugimoto model in [14, 15], also 
at zero temperature as in the field theory computation. 

In this letter we find supporting evidence for this QCD proposal by consid- 
ering a suitable model within gauge/gravity duality at finite temperature. The 
holographic model considered is ST/ (2) Einstein- Yang-Mills theory in (4 + 1)- 
dimensions. In the presence of a background B field, we do indeed find an 
instability, corresponding to the formation of vector hair in the gravity the- 
ory. In addition to the motivation from QCD described above, the result is 
also of more general relevance in the context of gravity and gauge/gravity du- 
ality. Whereas the QCD result has been obtained at zero temperature, here we 
consider the finite temperature case. 

In fact, the model we consider is a simple modification of the p-wave holo- 
graphic superconductor, first discussed in [5], in which there is a vector conden- 
sate in the background of a chemical potential and density turned on by con- 
sidering a non-trivial radial profile for the temporal component A t of the gauge 
field in the gravity theory. In addition, a p-wave superconductor for which the 
dual field theory is explicitly known has been constructed in [16, 17, 18] by 
embedding a probe of two coincident D7-branes, corresponding to two flavours, 
into the AdS-Schwarzschild black hole background. This string-theoretical top- 
down approach lets one identify the SU(2) chemical potential as an isospin 
chemical potential and the condensate as a p meson. A p-wave holographic 
model which includes the back-reaction of the gauge field on the metric in the 
SU(2) Einstein- Yang-Mills model has been considered in [19, 20]. 

Here we show that a similar condensation occurs also when a spatial com- 
ponent A y of the background SU(2) gauge field, corresponding to an external 
B field, has a non-trivial profile instead of the temporal component A t in the 
isospin chemical potential case. We find the corresponding holographic phase 
transition at finite temperature. For this we analyse the fluctuation quasinor- 
mal modes and find that they display an instability. We confirm this result 
by performing a Schrodinger potential analysis. Moreover we comment on the 
expected form of the condensate, the explicit calculation of which we leave for 
the future. 



lr This effect is similar to the Niclscn-Olscn instability for gluon condensation [10] and to 
W-boson condensation induced by a strong magnetic field [11, 12, 13] 
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2. Holographic Setup 



We consider SU(2) Einstein- Yang-Mills theory in (4+l)-dimensional asymp- 
totically AdS space. The action is 



/< 



_L (fl _ 2A) __L CF - 



+ Sbdy , (1) 



where K5 is the five-dimensional gravitational constant, A = — is the cosmo- 
logical constant, with L being the AdS radius, and g is the Yang-Mills coupling 
constant. The SU{2) field strength F« v is 

F^ = d^-d v A; + e abc A^Al, (2) 

where /j,, v = {t,x,y,z 7 u}, with u being the AdS radial coordinate, and e abc is 
the totally antisymmetric tensor with e 123 = +1. The A® are the components 
of the matrix- valued gauge field, A — A^^Ax^, where the r a are the SU(2) 
generators, which are related to the Pauli matrices by r a = a a /2i. 5*bdy includes 
boundary terms that do not affect the equations of motion, namely the Gibbons- 
Hawking boundary term as well as counterterms required for the on-shell action 
to be finite. 

The Einstein and Yang-Mills equations derived from the above action are 



4 / 1 \ 



(3) 



where the Yang-Mills stress-energy tensor T M „ is 



T in , = ^(f^F^-^F^F"" ) . (r,) 



In order to simplify the calculation we will from now on exclusively consider the 
probe limit K5/5 — > 0. Thus the metric at finite temperature is given by an AdS 
Schwarzschild black hole 

d.s 2 = ^ (-/( U )di 2 + dx 2 + dy 2 + dz 2 + , (6) 

with blackening factor f(u) = 1 — u A /u^. The Hawking temperature of this 
black hole is given by T = 1/nUh- Due to the probe limit we now only have to 
consider the SU(2) gauge fields living in this fixed background. 

In this paper we introduce a constant flavour magnetic field B = F^ y which 
solves the Yang-Mills equation of motion (4). Choosing the magnetic field in 
the z-direction breaks the rotational symmetry 5*0(3) down to £0(2) while 
choosing the third flavour direction breaks the SU(2) symmetry down to U(l) 
symmetry generated by r 3 rotations. We call this U(l) symmetry 1/(1)3. We 
may choose a gauge where only Ay = xB is non-zero (see Appendix A) . 
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3. Perturbations About the Equilibrium 

A system which is close to equilibrium can be described by linear response 
theory. There the effect of an external perturbation is given by the convolution 
of the retarded two-point function with the source of the perturbation. Due to 
the Cauchy integration formula, the response can be written as a sum over the 
poles of the two-point function (see e.g. [21]). These poles can be identified 
with the different quasinormal modes of the black hole [22]. 

Quasinormal modes of a black hole are distinct perturbations of the black 
hole solution. They can be understood roughly as resonances of the black hole. 
However since the energy of the perturbation can leak into the black hole, these 
fluctuations arc not normal modes and thus have been dubbed quasinormal. 
Each quasinormal mode's corresponding frequency consists of a real and an 
imaginary part. As for the damped oscillator, the real part of the frequency es- 
sentially determines the energy of the fluctuations, while the imaginary part is 
responsible for the damping. Writing the time dependence as e -1 "* we note that 
a relaxation towards equilibrium can only happen if all the quasinormal modes 
lie in the lower complex half plane. In AdS spacetimes the quasinormal modes 
satisfy the following boundary conditions: at the horizon they are purely ingo- 
ing, whereas at the conformal AdS boundary they have an asymptotic behaviour 
that corresponds to a normalisable mode. 

Let us now study perturbations of the gauge fields about the equilibrium 
given by the magnetic field A y . As usual we fix a gauge by setting a£ = 0, which 
however leaves some residual gauge freedom. It is useful to consider combina- 
tions of the fluctuations which do not transform under the residual gauge 
transformations. These combinations are the physical modes of the system. 
First we define E^ — a* ± ia^ which are in the fundamental representation of 
i7(l)3 while are of course in the adjoint (see e.g. [23]). Since fluctuations 

cannot interact with the magnetic field, their equations of motion are the 
same as for zero magnetic field. Thus we cannot observe any instability there 
and we will consistently set = 0. We can also consistently choose that the 
fluctuations have a dependence only on t, x and u. 

Now we combine the fields E in such a way that they are covariant under 
the residual gauge symmetries, i. e. they only transform as a fundamental field. 
Using the Fourier ansatz e -1 "*, this leads to the six gauge covariant fields, (see 
Appendix A) 

E± x = x 2 B 2 Et ± i {xB8 x Ef - BEy) ,E± 2 = ±xBEf + cjE^ , E$ = E$ . 

(7) 

The equations of motion for E^ x y , as derived from (4), are given by the 
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coupled set of equations 

=dlE+ \d u Ef { -^Ef + - f dlEt ^Ef + Jd x Ef , (8) 
=dlEi + duEi { -^Ef + j(l xd x )Ef 



1W ^J- 



f 2 d x Et + j,Ef , (9) 
=d 2 u E+ + (-1 + d u E+ j (2 + xd x )Ef + jd 2 x E+ 

+ ^ E t + j2 E i' ( io ) 

= - \Bxfd u E+ - fd x 8 u E+ - iu)d u E+ , (11) 

where the equations of motion for E^~ x ^ y are obtained by complex conjugation. 
Notice that it is consistent to switch off Ef~ since (8) then coincides with the 
constraint (11) and the gauge covariant fields E^ 2 disappear. The equations of 
motion for Ef decouple from the above system since they are the only transverse 
modes and thus can be studied on their own. 



4. Finding the Instability 

In order to find the instabilities of the system we determine the spectrum 
of the quasinormal modes and determine the critical value of the magnetic field 
B for which the imaginary part of the quasinormal modes becomes positive. 
Usually the crossing to the upper complex half plane occurs at the origin (see 
e. g. [5, 24]) which we also expect here. By considering the equations of motion 
(8)— (11), we observe that Ef decouple if u> = 0. Thus we can look for normal- 
isable modes of these fluctuations for cu — separately. Since we cannot find a 
normalisable mode for any magnetic field, we can conclude that the fluctuations 
Ef- do not trigger any instability. Therefore we consider from now on Ef- = 
which is consistent. 

Also, the equations of motion for Ef decouple from all the other equations 
of motion since they are the transverse modes. We numerically checked that 
these fluctuations do not introduce an instability either. Thus an instability 
can only come from the E x y components, the quasinormal modes of which are 
studied in the following. 

4-1. Analysing the Quasinormal Modes 

For Ef = it can be shown that the set of equations (8)— (11) reduces to 
one equation of motion for Ef, 

= dlEf + ^dlEf duEf ^d x Ef + - Ef , (12) 
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while E+ is given by the constraint 

K = ^ 9 * E 2- ( 13 ) 

Usually it is favourable to write down the equations of motion in terms of the 
gauge covariant fields (7). That is also possible in this case, but here we obtain 
the solution for the gauge covariant field in terms of the solution for fi+ by 
using the constraint (13), 

El, = ~(dl- \d x x 2 B 2 ^ Et . (14) 

From this it is clear that the dependence on the radial coordinate u is the same 
for both fields, and only the dependence on x may change. 

The solution for can be found by the separation of variables E+ = 
X{x)U{u) and its equation of motion (12) becomes 

= X" - -X' - (C + B 2 x 2 ) X , 

where C is the separation constant. Using (14), the equation for X and the fact 
that E^ is separable, we find that 

Eli = -CEt ■ (16) 

This shows that the physical fields E^ 1 and E^ have the same X- dependence 
and thus are equivalent up to a constant. In solving for X, we will have to impose 
boundary conditions that ensure that the vanish asymptotically and are 
regular. As we will find out, they have an overall factor of exp(— Bx 2 /2). 

The equation for X can be solved analytically. After the definition X = 
exp(— Bx 2 /2)Y and y — Bx 2 we obtain Rummer's equation 

yY"(y) + (b-y)Y'(y)-aY(y)=0, (17) 

with b = —1/2 and a = (C — B)/AB. The solution is given by the linear 
combination of the confluent hypergeometric function of first kind M(a,b,y) 
and of the second kind U(a, b, y). The function U can also be written in terms 
of A4 such that the complete solution is given by the linear combination 

Y = dM(a, b, y) + C 2 y 1 - b M(l + a-b,2-b,y). (18) 

In order to have a vanishing solution for E^ 1 (and thus for X) at x — > oo, Y is 
only allowed to grow slower than the factor of exp(y/2) since X = cxp(— y/2)Y. 
Looking at the asymptotic expansion for large y, 



6 



we see that Y does grow too fast unless T(a) diverges, which is to say that 
a = —m for a non- negative integer to. Indeed, when a = —m then M.(a, b, y) 
reduces to a Laguerre polynomial Lm\y), which has m positive roots. We 
expect that the free energy of the system increases with the number of the 
roots. Since we are looking for the lowest energy solution and the second solution 
in (18) always has a root at y = 0, we obtain C2 = 0. The first solution has no 
zero if a = 0. This determines 

C = B and X = -e~ Bx2 ' 2 , (20) 

where the normalisation is set by X(0) = —1. We thus find that 

= Bc- Bx2 / 2 U(u) . (21) 

Unfortunately the equation for U can only be solved numerically. Using a 
shooting method we may find its quasinormal modes and determine the depen- 
dence of their frequencies on the magnetic field. In order to have dimcnsionless 
quantities we define 2$ = Bu\ = B/(irT) 2 and to = uju h = uj/itT. For *B = 
the quasinormal frequencies are given by lt)„ = n(±2 — 2i) [25]. The dependence 
of the quasinormal frequencies on the magnetic field is given in figure 1. We 
observe that for the magnetic field *8 ~ 5.15, the quasinormal frequency moves 
into the upper half plane, which signals an instability. 

4-2. Tachyonic Mode 

There are two further ways of qualitatively understanding the formation of 
an instability that leads to a superconducting condensate. The first is by looking 
at the form that the on-shell action takes in our setup. To see an instability 
we actually only need to switch on the fields Ay = xB, a\ and a? y . To second 
order, the relevant part of the action becomes 

S(2) = ~ f — {2,/ V 1 (dtal) 2 + 2g u gyy (d t al) 2 + 2g uu g™ {d u a} x f 



+ 2g^ g yy {d u a 2 y ) 2 + 2 9 xx g yy [B 2 + {d x a 2 y f 

+ 2Ba\a 2 y - 2Bxa\d x a 2 y + B 2 x 2 (a*) 2 ] j . 



(22) 



Solving the x- dependence, we obtain _E+ = e ^ Bx V(t,u) and equation (13) 
gives £7+ = —iE+ which translates to a x = —a 2 = c~2 Bx Re V. This means 
that, to second order, the on-shell action is 



S {2) = y dtdydzdu 



2g tt (gxx + g yy) (d tRe yf 



+ 2g uu (g xx + g y y) (a u Rc Vf - ABg xx g y y (Re V) 2 



(23) 
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Figure 1: The movement of the quasinormal frequencies at increasing f B. The two lowest 
order frequencies start at tt) = ±2 — 2i. They move towards the imaginary axis, merge, then 
split again. The merging happens when <8 ~ 1.92. One pole then moves downward, away 
from and outside the range of the figure. The other, which is drawn, moves towards and 
eventually crosses the real axis when 2$ = 23 c fa 5.15. 

where we have integrated out the x component, absorbed a factor of (n/B) 4 into 
V for the canonical normalisation and dropped a constant term. The important 
thing to notice is that the gauge field has picked up an effective mass term, 
m 2 ff = —4Bg xx g vv , and it becomes more and more negative as B increases. This 
suggests a tachyonic mode leading to the spontaneous formation of a condensate. 
In contrast to the instability at finite chemical potential [5] , the instability here 
does not arise due to the presence of a horizon, i. e. g is not involved in the 
negative effective mass term. 

4-3. Schrddinger Potential Analysis 

To understand the instability better, we can rewrite the equation of motion 
(12) in the form of a Schrddinger equation (see e.g. [26, 27]), 



where E — to 2 and ip oc E£. The tortoise coordinate R and the Schrddinger 
potential are given by 



R{u) cannot be inverted analytically and so we have to study V s (R) numerically. 
The profile of the Schrddinger potential for different values of the magnetic field 



d 2 R t[> + V s iP = Et{;, 



(24) 



R = — (arctanh u + arctan u) and V s 



(l-u 4 )(3-4Bu 2 + 5w 4 ) 
4u 2 



(25) 
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Figure 2: The Schrodinger potential at different values of OS. It is helpful to compare this with 
figure 1. From top to bottom, the curves show the potential at the values 23 = 0, 1.92, 5.15 
and 10. The first curve shows the potential when there is no magnetic field, the second one is 
for the value of 23 at which the poles merge on the imaginary axis, the third is for the value 
at which the first pole crosses the real axis, and the final value is to see what the potential 
does at higher 23. We see that a potential well forms with negative potential energy. The 
value 23 = 1.92 is roughly where the curve first dips below zero, which is where the poles in 
figure 1 merge. 

03 is given in figure 2. For all values of the magnetic field we observe an infinite 
wall at R = due to the boundary of the AdS space and a vanishing potential at 
R — > oo due to the black hole horizon. In between a potential well opens up as 
the magnetic field is increased. For a large enough magnetic field, this potential 
allows for bound states with negative energy. The first such bound state implies 
that a quasinormal frequency has crossed the real line and developed a positive 
imaginary part. This bound state therefore first forms at the critical magnetic 
field that produces an instability. 

The reason why a bound state with negative energy implies an instability is 
due to the ingoing boundary condition at the horizon. This boundary condition 
implies that the time- independent part of the mode goes like e lujR . Outside 
the potential well, the mode has to decay exponentially towards the horizon. 
Combining these two constraints means that u> must have a positive imaginary 
part. Thus the appearance of the bound state in the Schrodinger picture agrees 
with the quasinormal frequency crossing the real axis. 

Using the Schrodinger potential we can also estimate the critical B-field Q5 C , 
above which the system gets unstable. Using a WKB approximation, we can 
estimate that a (zero-energy) bound state will appear for 

y* oo p 1 j 

tt= / dR^V^R)= jv^W, (26) 

where n £ N. Rq is the zero of V S (R) and thus V S (R) is negative for R > Rq. 
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The corresponding zero of V s as a function of u is called uq, 



u = ^- |V4«8 2 - 15. (27) 

Note that V s (u) and therefore also wo depend on 03. The minimal value for 03 
where the potential develops a zero is at 03 = -^p ~ 1.94. This coincides with 
the value at which the quasinormal modes merge at the negative imaginary axis 
up to an error of 2%. 

In order to obtain the critical magnetic field we need to know when a bound 
state can be formed. We are only interested in the first bound state and therefor 
set n = 1. Solving the integral equation (26) numerically for 03 we obtain 03 c s=s 
5.3514. We see that the analysis of the Schrodinger potential qualitatively agrees 
with the precise study of the quasinormal modes above up to an error of 4%. 



5. Conclusion 

In this letter we have found an interesting instability of black holes in asymp- 
totically AdS spacetimes. In particular we consider an Einstein- Yang-Mills sys- 
tem with gauge group SU(2) and with a negative cosmological constant such 
that AdS is a solution. If we neglect the back-reaction of the gauge fields on 
the background metric, a solution of the equations of motion is given by an 
asymptotically Schwarzschild black hole. 

In addition wc turn on a magnetic field F xy — Bt 3 , where t 3 is one of 
the three generators of SU(2). By investigating the fluctuations of the gauge 
field around the Schwarzschild black hole we find that the system gets unstable 
above a critical value 03 c w 5.15 of the magnetic field. In particular we identify 
the quasinormal mode which crosses the real axis into the upper half plane. 
This result is confirmed in two different ways. First, we calculate an effective 
action for the fluctuations and show that these fluctuations get a negative mass 
squared, which drops below the Breitenlohner-Freedman bound for large enough 
magnetic fields. Second, we rewrite the equations of motion into a Schrodinger- 
like wave equation with an effective potential V s and show that V s has a bound 
state for large enough magnetic fields. The bound state signals the instability 
of the system. 

Note that here we have turned on a magnetic field for one of the three 
generators of SU(2) and not for a U(l) gauge group. One might think that this 
is not so common in gauge/gravity duality. However, as we now discuss, a similar 
structure [14, 15] also arises in the Sakai-Sugimoto model [28]. To see this, let 
us consider the Sakai-Sugimoto setup with two flavours (with electric charges q\ 
and (72). The flavour degrees of freedom are described by two pairs of D8/D8- 
branes. Part of the dynamics of the flavour degrees of freedom is therefore given 
by a non-Abelian gauge theory with gauge group U(2)l x U (2)r. Note that the 
gauge group on the gravity side corresponds to the dual flavour symmetry on the 
field theory side. In the Sakai-Sugimoto model the electromagnetic gauge group 
with generator Q = diag(gi, g 2 ) is realised as a subgroup of U(2) l*U {2)r. Thus 
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in the case qi ^ q?, the magnetic field B em has a baryon and isospin component, 
i.e. 

cm t 3 - (28) 

Although we aim to study the effect of the U(l) magnetic field B em , we have to 
consider a non-vanishing magnetic field for the isospin component. Note that 
already in [14, 15], evidence was given that the vacuum at zero temperature is 
unstable for large enough magnetic fields. The results presented here indicate 
that such an instability arises not only in the Sakai-Sugimoto model but also in 
the simplest possible toy model at finite temperature for which we can explicitly 
study the quasinormal modes. 

The instability discussed in this paper gives rise to interesting questions on 
which we plan to work in the future. What is the true ground state of the 
system and what are its properties? Is it a superfluid state? Can we realise this 
ad hoc toy model in string theory? By addressing these questions, we expect 
in particular to obtain further information about the QCD vacuum instability 
at finite B field discussed in [6, 7, 8, 9]. As our results suggest, a more general 
form of this instability may occur also at finite temperature in non-confining 
field theories. On the gravity side, we expect the condensate to be a suitable 
combination of the gauge field components A\ and A y . 

It will also be interesting to repeat the analysis of this paper in four di- 
mensions on the gravity side. This will allow for a comparison with [29], 
where instabilities of general spherically symmetric solutions of four-dimensional 
Einstein- Yang-Mills-Dilaton theories with a constant magnetic field strength on 
the sphere have been found. Moreover, this will also allow for a comparison with 
[30], which raises the interesting question whether the instability found here is 
related to the existence of magnetic monopoles. 
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Appendix A. Constructing the Gauge Covariant Fields 

In this appendix we show how to fix the gauge of the SU(2) field in such a 
way that only a 17(1)3 gauge freedom remains. We then find the fluctuations 
that transform covariantly under the remaining Z7(l)3- 

For the background field A° without fluctuations, the gauge transformation 
e* A ( x ") acts infinitesimally on A a as 

5AI = df,A a - e a6c A 6 A£ . (A.l) 
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As usual we fix the gauge such that A u = 0. This leaves A a (t,x,y, z) as the 
residual gauge symmetry. We choose to have the magnetic field B = F 3 y . 
From (2) we see that this means we can fix A a x = 0, Ay 2 = and Ay = xB. 
These choices give that F" = d x Ay — B5 a3 . The choice that A% — reduces the 
remaining gauge symmetry to A a (t, y, z), and the conditions on A® fix A 1 ' 2 = 
such that only A 3 (t, z) is left. The SU{2) has been broken to C/(l) 3 . 

Now we consider the fluctuations a^(i, x, y, z, u) about the background field 
Ar\. The full field A® + a° transforms under the infinitesimal gauge transfor- 
mation A a + A a . We treat A a as an infinitesimal perturbation of A a . Then the 
fluctuations af. transform as 

5al - d^\ a - e abc (A b a^ + \ b A^) , (A.2) 

while the transformation for A^ remains the same as before. Again we fix a u = 
which leaves X(t, x, y, z) as the residual gauge symmetry for the fluctuations. We 
define the fields E ± — a 1 ± ia 2 . The gauge transformations of the fields E ± and 
a 3 to first order are 

5E± = d y X ± T iA 3 S± ± iA 3 y X ± , (A.3) 
Sal = d ^ > 

where A ± = A 1 ± iA 2 . Notice that the fields Err transform in the fundamental 
of the U(l) 3 . We now consider only a dependence of t and x and use the 
Fourier ansatz e~ lwt . The fields transforming covariantly under the U(l)s gauge 
symmetry are then 

Ef A = x 2 B 2 E~f ± i (xBd x E$ - BEy) , E 3 = d x a\ + ]ua 3 x , 

E± 2 = ±xBEf + urE± , E 3 A = a 3 y , (AA) 

E T = E t i E T,2 = a l ■ 

Since the B field does not interact with fields in flavour direction 3, these U(l)s 
covariant fields are trivial. However the fields in direction 1 and 2 couple and 
take a non-trivial form. 
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